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Abstract 

We present a simple derivation of the matrix Riccati equations governing 
the reduced dynamics as one descends from the group U (N) describing the 
Schrodinger evolution of an iV-level quantum system to the various coset spaces, 
Grassmanian manifolds, associated with it. The special case pertaining to the 
geometric phase in iV-level systems is described in detail. Further, we show 
how the matrix Riccati equation thus obtained can be reformulated as an equa- 
tion describing Hamiltonian evolution in a classical phase space and establish 
correspondences between the two descriptions. 
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1 Introduction. 



There is currently considerable interest in studying various properties of quan- 
tum systems with finite-dimensional state spaces, for instance in understanding 
entanglement in the context of Quantum Information and Quantum Computa- 
tion pQ. Perhaps surprisingly, there seem to be many things yet to be learnt 
in such systems. From the standpoint of fundamental Quantum Mechanics as 
well, there has been longstanding interest in for example extending the concept 
of Wigner distribution to finite dimensions and its role in quantum state esti- 
mation [5] . In this and other contexts one finds that many new features specific 
to the number of dimensions appear[3]-[5| that often have no counterpart in the 
infinite dimensional case . 

For two-level systems, N = 2, it is well known that the Poincare-Bloch 
sphere § 2 provides an excellent realization and practical tool for dealing with 
pure (as well as mixed) quantum states [5] . Some efforts to extend the Poincare 
sphere construction to N > 3 also exist [7]. In another direction, there has 
been a systematic programme of 'unitary integration' methods [SJ, which also 
generalize the Poincare sphere concept to higher dimensions, and lead to meth- 
ods of projecting quantum dynamics at the iV-level Hilbert space governed by 
a Schrodinger equation to various 'base spaces' of lower dimension. From a 
mathematical point of view, the central ideas and structures underlying these 
'reductions' are similar to those known and extensively developed, over a long 
period of time, in the context of classical dynamical systems with the objective 
of seeking dynamical superpositions in non linear evolution equations [9]. In a 
recent work of Uskov and Rau [10] , one begins with the group SU (N) acting on 
the state space of an TV-level system, the Schrodinger unitary evolution opera- 
tor being also an element of SU (N). For any partition: N = n\ + n 2 , one has 
the subgroup SU (m) x SU (n 2 ) and the coset space SU (N) /SU (m) x SU (n 2 ), 
which is a Grassmanian manifold and which functions as the base manifold of 
a fiber bundle, the total space being SU (N). Using matrices outside of SU (N) 
generated by nilpotent matrices in a well-chosen manner, a parametrization of 
the base manifold by a set of n\ ■ ri2 independent complex coordinates is set 
up. It is then shown that the original Schrodinger evolution projects down to a 
system of matrix Riccati equations for these base-space coordinates, and some 
connections to geometric phases [IT] are indicated. 

The aim of the present work is to give a treatment of this problem of re- 
duction of Schrodinger evolution to various base spaces in a manner that works 
throughout within the unitary group U (JV) intrinsic and natural to Quantum 
Mechanics, and to obtain the Riccati equations in a rather elementary manner. 
Thus the use of 'nilpotent generators' is entirely avoided, and in addition the 
region of U (N) covered by the usual complex variables is easily seen. Follow- 
ing this, in the case n\ = N — 1, ni = 1 appropriate for discussing evolution 
of pure states by the Schrodinger equation, it is shown that the base space can 
be conveniently viewed as a classical phase space, and the Riccati equations 
expressing quantum dynamics then appear as classical Hamiltonian equations 
with a suitable Hamiltonian function. This framework is then used to describe 
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in detail the structure of pure state geometric phases for the original TV-level 
system: one sees the extent to which such phases can be expressed in purely 
classical terms, and also particularly clearly that they are ray-space quantities. 

A brief summary of the present work is as follows: In Section 2 we derive 
the Riccati equations associated with various coset spaces of the group U(iV) 
governing the Schrodinger evolution of an N level quantum system without 
resorting to nilpotent operators. In Section 3 we apply the general formalism 
of Section 2 to a special case appropriate to the geometric phase. Here we also 
show how the resulting equations can be reformulated in purely classical terms 
and compare the two perspectives on the geometric phase. Section 4 contains 
our concluding remarks. 



2 TV-Level System Dynamics, Coset Spaces, Ric- 
cati Equations. 

Consider an iV-level quantum system with a Hermitian Hamiltonian matrix 
H (t) which may be time-dependent. The unitary evolution operator U (t) is an 
element of G = U (N) obeying the equation: 

iU(t) = H(t)U(t), U{t )=I. (1) 

When H is time- independent, one sees easily that the solutions of this 
equation are given by the one-parameter group generated by iH acting on G 
from the left. 

The right action of G on itself acts transitively on the family of solutions. 
Therefore, all the solutions are simply obtained by acting on the subgroup gen- 
erated by iH. This circumstance gives rise to the nonlinear superposition rule 
present in Lie-Scheffers systems [9]. 

As discussed below, when we consider 'decompositions' of G, we may relate 
our equations of motion with other evolution equations on the space defined by 
the decomposition, the "reduced space". On this reduced space the evolution 
equations are in general nonlinear and acquire the form of Riccati- type equation. 
Many 'reduction procedures' are available in the literature, see for instance |12j 
for some applications. What we consider here are specific instances applied to 
relevant physical situations which arise in quantum mechanics. 

To 'decompose' the unitary group G, let us consider any partition: N = 
ni + ri2 and the subgroup: 

H = U (m) x U (n 2 ) C G, (2) 

with the factor U (m) acting on the first ni dimensions, U (712) on the rest. The 
idea is to describe the coset space G/M 'nicely' and obtain from Eq. {1} by 
'projection' an evolution equation on it. Write a general matrix in G in block 
form as: 

A B 
C D 



U 



(3) 
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where A,B,C,D are ri\ x x n 2 ,n 2 X n\,n 2 X ^2 dimensional matrices 

respectively. The condition Wll = 1 becomes: 



A^A + tfC = 1, A t B + C t £) = 0, D^D + B^B = 1. 
Under right multiplication by an element 



Ui 
U 2 



in H, the changes in A, B, C, D are: 

A — » At/i, B — > -BC/2, C^CUu D^DU 2 



(4) 



(5) 



So, in case A and -D are non-singular, C A~ x and BD~ X are invariant, i.e. they 
are constants over each coset, and so are CA> and SD^ 

Let us now limit ourselves to the subset of G in which A, D are both non- 
singular. Then one can get a nice coset representative as follows: using the 
freedom of right multiplication by elements in H, we can use the polar decom- 
position of A and D to bring them both to Hermitian positive-definite forms: 



coset representative = Uq = 



A B 
Co D 



A\ = Ao> 0, D\ = Do > 0. (6) 



The coset invariant mentioned above can then be defined as an m x n 2 complex 
rectangular matrix: 

Z = B D \ (7) 
We can now use Eq.(|3]) to express Aq, Bq, Co, Dq in terms of Z: 



Ao = T- 1/2 , Bo = ZT~ 1/2 , 



c a = -z^T 1 1/2 , D = r 2 



1/2 



(8) 



where T 1 = I + ZZ^ and T 2 = I + Z^Z are n\ x n\ and n 2 x n 2 positive 
hermitian matrices intertwined with each other through Z: Y\Z = ZT 2 . Thus, 
in the region of G defined above, a general matrix is U in Eq. can be written 
as: 



U 



A Bo 
Co Do 



Ui 

u 2 



(9) 



giving A = AqU-l, B = B U 2 , C = C Ui, D = D Q U 2 . Using this in the 
evolution equation ([1]) and writing H (t) in block form: 



we get: 



A Bo 
Co D 



H(t) 



Ui 

u 2 



ffiW V{t) 
V(t? H 2 (t) 



Hi V 
Ut H 2 



(10) 



A B 
Co Do 



Ui 

u 2 



4 



i.e.: 



i 


A 


B 






V 




A Q 


Bo 


— i 


A 


Bo 




UxU^ 1 







Co 


Do 






H 2 




Co 


Do 




Co 


Do 







W 1 



which yields 



iA = HxAo + VCo-iAoU^ 1 , 
iBo = H 1 B Q + VDo-iBoU 2 U^ 1 , 
iC Q = V^A + ^Co-iCoU^ 1 , (11) 
iD a = V^Bo + HsDo-iDaUiUi 1 . 
Using the second and fourth of these we get an 'autonomous' equation for Z: 

iZ = iBoDa 1 - iBoD^DD^ 1 

= HiZ + V- iBoU 2 Uz l D a - 1 - Z \V^Z + H 2 - %DoU 2 U^ l D^ 

= V + HiZ - ZH 2 - ZV ] Z , (12) 

where, for brevity, the ^dependencies of V, Hi,H 2 , have been omitted. 

This is a matrix Riccati equation for evolution on the coset space. The com- 
plete quantum evolution ([T]) involves also equations for U± and U 2 . This deriva- 
tion seems simpler and more direct than others [10j . and shows the appearance 
of the Riccati structure rather clearly, including the result of the projection. 

3 Connections to the Geometric Phase. 

To connect up to the Geometric Phase, we limit to the choices n\ = N — 1 , n 2 = 
1. So in the Hamiltonian matrix H (t) of ([ID]). Hx is (N - 1) x (N - 1), V is an 
(N — l)-component column vector, and H 2 is a single real quantity. We recall 
here some notations: 

• HS N ^: complex iV-dimensional Hilbert space consisting of vectors ip,ip'. 

• Bn : unit sphere in TL^ of real dimension (2N — 1) identified as a coset 
space U (N) /U (N - 1) = G/H , H = U(N-1). 

• IZ n :space of unit rays, of real dimension 2 (N — 1) identified as a quotient 
and coset space B N /U(1) ~ G/H — U (N) /U (N - 1) x U (1) ; H = H x 
U{1) = U{N-1) x U{1). 

We now consider four aspects: 
(a) (Local) coordinates over TZn ■ 
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Let us limit ourselves to that part of T-6 n \Bn in which the last, N — th 
component of ip € 7i^ N ^ is non-zero. Then, for a vector ip € Bn and its image 
p € TZn we can say: 



yl/2 > 



, z € W^" 1 ), 7 = 1 + zh; ip = e za ip (2) , 



P = tp^ = "00 (z) "00 (z) f , 1P0 (z) 



(13) 



So z e TiS N l > becomes a system of coordinates for the base space TZn = 
B N /U (1) ~ G/H = U (N) /U (N-l)x U (1). When useful we will later write: 

z r = q r + ip r , r = l,2,...,N-l, (14) 

so q r ,p r are 2 (TV — 1) real independent local coordinates on ray space. 

(b) Schrodinger equation on TC^ N ' to Riccati equation on TZm 
Consider the Schrodinger equation for a vector ip (t) S TlS N ^: 



iip (t) = H(t)i>(t) = 



B x V 



1>(t) 



(15) 



Separating the upper (N — 1) components of ip denoted by £ from the N — th 
one 77: 



(16) 



Eq. p5p may be rewritten as: 

i£ = H 1 € + Vri, iri = V^ + H 2 n. (17) 
If we now set: z — £/r/, in analogy to Eq.|7]), we get its evolution equation: 



iz = i^/rj — i^rj/rj 2 

= V + H x z- zH 2 - zV f Z. 



(18) 



This is again a Riccati equation, like (|12[) . but now z is an (N — l)-component 
column vector, as is V, and H 2 is a single real variable. Now we have obtained 
the Riccati structure for a single solution of the Schrodinger equation, not using 
the entire unitary evolution operator. So the quantum-mechanical ray space 
evolution is given by the non- linear Riccati equation (fT5|) . 

(c) Reformulation in classical phase space form. 



We now show that Eq. (fl8|) can be written in completely classical form ouTZn 
regarded as a phase space. We start from Eq. (fT3"|) and define first a one- form # 
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and then a symplectic two- form lo on 1Z n as follows: 

9 Q = -iip (z) f d^ (z) 



and: 



where 



= — Im z^dz 
7 

= — (q r dp r — p r dq r 
7 



-1/2 







+ 



-1/2 



wo = dd = —idipo (z)^ A dipo (z) 

2 2 
= -dq r Adp r (q r dq r + p r dp r ) A (q s dp s - p s dq s 

7 7 2 



2 I dq r dp r I A 



— M rs L r; 



dq s 
dp s 



L rs = -L sr = (q r p s - q s p r ) /2j , 



M rs = M s . 



27 



(q r q s +PrPs) /27 2 . 



The two-form ujq is indeed non-singular, and the inverse of the 2 (N 
2 (N — 1) matrix above can be computed: 



where 



L 


M 


-l 


X 


Y 


M 


L 




-Y 


X 



X rs = -X sr = 27 (q r p s - q s p r ) , 

Y rs = Y sr = -27 (S rs + q r q s + p r p s ) , 

and we can read off the fundamental PB's on ray space in q — p and z — z 

Ur,q S } = {Pr,Ps} = l(qrPs - qsPr) /2 = X rs /4:, 

{q r ,p s } = -7 (S rs + q r q s +PrPs) /2 = Yrs/4, 

{z r ,Z*} = ij (S rs + Z r Z*) , {z r , Z s } = {z* , Z* g } = 0. 

From here the following useful PB's can be obtained: 



{z r ,j} = i-f 2 Z r , {Zr,7 1 } 



{Zr,Z* s h} 



lS r 



{Zr,f(z)/l} = -iZrf(z). 
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Here / (z) is analytic in the z r 's. Using all this, we can show that the quantum- 
mechanical Riccati evolution equation (|18[) can be written as a purely classical 
Hamiltonian evolution on TZn with the P_B's given in 



-i (V + H x z - zH 2 
{z 7 H(z,z^)}, 



(26) 



where 



= -(H 2 + z^V + V i z + z^H 1 z)/- f 
1 



-I | 2 t 

7 



i?i V 
V* H 2 



-Vo (z)^H^ (z) 



(27) 



The motion of z in TZ n, induced by unitary Schrodinger evolution in TL^ N \ 
is thus a continuous classical canonical transformation. 

(d) Geometric phase and its 'classical' aspects 

We consider Schrodinger evolution as given by Eq. (|15[) . with no appeal to 
the adiabatic approximation. The Riccati equation (|18[) for z, as we have seen, 
is a consequence of Eq. p5[) . and from it we find for 7: 

1 {t)=i 1 {t)(y^z{t)-z{t)W), (28) 

where in general V — V (t). Comparing Eqns. (fT5) ) and p^|) gives: 

e = e i VV7, ^e M /V7- (29) 

Using these and (|28[) in the r\ equation of motion in (|17|) gives the evolution 
equation for the phase a (t) : 



a = -H 2 - i (Ft z + z iv) 



(30) 



( Both Eqs. (|2"8"l) ( j3"0"l) appear in [TU]). In summary, the Schrodinger equation (TT5j) 
for -0 (t) amounts to an autonomous Riccati equation (jlSp for z, re-expressed in 
(|25)) in classical phase space form and equation of motion ([30]) for a (t) where 
the right-hand side is a independent. 

If now tp (t) is any (not necessarily cyclic) solution of the Schrodinger equa- 
tion (|15|) between the given times t\ and t 2l from general theory it is known 
that the geometric phase is the difference of two terms, a total and a dynamical 
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phase: 



^Pgeorn ^Ptot ^Pdyn j 



f tot = arg (-0 (ti) , ip (t 2 )) 

= a (t a ) - a (t x ) + arg (l + z f (t x ) z (t 2 )) , 




ti 

*2 



dttpo (z (t)y H (t) (t) 



ta 



dtH{z{t),z{t)*) . 



(31) 



ti 



The first part is calculable using ([30]) and depends only on ray space quantities; 
the second part is the time integral of the classical Hamiltonian along the ray 
space trajectory z{t). Thus the geometric phase is seen to involve only ray 
space quantities, as it must, and is expressed as far as possible in terms of the 
classical Hamiltonian Ti (z, z*) via the dynamical phase. In special case of cyclic 
evolution, z (t\) = z (tz), the total phase simplifies to just a (£2) — en (ii). 

If we adopt the kinematic approach in which there is no use of the Schrodinger 
equation, the results our similar except for a shifting of terms. We now consider 
directly a (smooth) parametrized curve C in Bm with image C in TZn'- 



C = 



C = {^j (s) = e^Vo (2 («)) |*i < a < s 2 } C B N: 
= {p (a) = ip (z («)) Vo (2 (s)) f |si < s < s 2 } C TZ N . 



(32) 



The latter is essentially a curve in the space of z: 



C = {z(s) \si < s < s 2 }. 



(33) 
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Then the geometric phase is again the difference is as usual: 

tfigeom [C] = (fitot [C] - (fidyn [C] , 

iptot[C] = arg (<ip (si) , ip (s 2 )) 

= a (s 2 ) - a (si) + arg + z (si) 1 z (s 2 )) , 

f d 

^djm [C] = Im / dsV> (s) f — ^ (s) 

= Im / ds^(s) t |ia(s)i/;(s) + e iQ(s) -^Vo(z(s)) 



a(s 2 ) - a(s x ) +Im / V^o (z (s)) t d^o (z (s)) 

si 

along C 

a(s a )-a(*i)+ / Q . (34) 



Here we used Eq. (|19l) . Therefore again we have a ray space quantity: 

Wgeom [C] = arg (l + z {s^ z (s 2 )) - J O . (35) 

c 

If C is a closed loop, then z (si) = z (s 2 ), the geometric phase is a purely classical 
"symplectic area": 



dC = : </? geom [C] = - / O = - / / wo, (36) 
where 5 is any two-surface with dS = C . 



4 Concluding Remarks. 

In the present work, we have attempted has to highlight the connections between 
the following in the quantum evolution of TV-level systems: 

1. Appearance of Riccati equations in as direct a way as possible to describe 
'unitary' evolution in the coset space G/U (ni) x U (n 2 ) where G = U (N) 
and n\ + n 2 = N; 

2. In the case n\ = N — l,n 2 = 1: to recast the Schrodinger equation for tp 
in H^ N ' as a Riccati equation on the ray space TZn plus an equation (|30| 
for the overall phase a 'driven' by this ray-space evolution; to express the 
ray-space evolution as a classical Hamiltonian evolution with PB's (|24l) 
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and Hamiltonian function H. (z); to study the structure of the GP in both 
Schrodinger evolution and kinematic case, bring in the classical symplectic 
structure as far as possible, and show that the GP is always a ray-space 
quantity. 

We hope the approach to Riccati equations developed in the present work, 
owing to simplicity and directness, will find useful applications in problems 
involving 'reductions' as illustrated here with geometric phases as an example. 
It will be interesting to consider the situation in which G is the unitary group 
associated with a composite Hilbert space of dimension TV = nm, with the 
Hilbert spaces of the subsystems being of dimension n and m respectively. The 
subgroup H will be the tensor product of U (n) and U (m), the so-called group 
of local transformations. In this way the reduced space will be related to the 
set of entangled states of the composite system and this may shed new light on 
the classification problem of separable versus entangled states. 
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